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We show the emergence of a strongly interacting Bose-Fermi mixture from a two-component Fermi 
mixture with population imbalance. By analyzing m situ density profiles of ®Li atoms in the BCS- 
BEC crossover regime we identify a critical interaction strength, beyond which all minority atoms 
pair up with majority atoms, and form a Bose condensate. This is the regime where the system 
can be effectively described as a boson-fermion mixture. We determine the dimer-fermion and 
dimer-dimer scattering lengths and beyond-mean-field contributions. Our study realizes a Gedanken 
experiment of bosons immersed in a Fermi sea of one of their constituents, revealing the composite 
nature of the bosons. 
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Fermions are the fundamental building blocks of mat- 
ter, whereas bosons emerge as composite particles. The 
simplest physical system to study the emergence of 
bosonic behavior is a two-component fermion mixture, 
where the composite boson is a dimer of the two different 
fermions. A dramatic manifestation of bosonic behav- 
ior is Bose-Einstein condensation, representing the low- 
temperature phase of a gas of bosons. One way to reveal 
the composite nature of the bosons is to immerse such a 
Bose-Einstein condensate (BEC) into a Fermi sea of one 
of its constituents. The degeneracy pressure due to the 
Pauli exclusion principle affects the structure of the com- 
posite boson, resulting in a zero-temperature quantum 
phase transition to a normal state where Bose-Einstein 
condensation is quenched. 

In this paper we observe this transition experimentally. 
We identify the regimes where a two-component Fermi 
gas can be described as binary mixture of bosons and 
fermions, and where the composite nature of the boson 
becomes essential. The validity of a Bose-Fermi (BF) 
description requires that all minority fermions become 
bound as bosons and form a BEC. We determine the 
critical value of l/kpi^a for the onset of superfluid be- 
havior in the limit of large population imbalance. Here 
a is the fermion- fermion scattering length and kpij is the 
Fermi wave number characterizing the depth of the ma- 
jority Fermi sea. Of course, for an equal mixture, the 
zero-temperature ground state is always a superfluid in 
the BEC-BCS crossover. It has been shown previously 
that a crossover superfluid can be quenched by popula- 
tion imbalance, also called the Chandrasekhar-Clogston 
(CC) limit of superfluidity [i,^. In this work we deter- 
mine the critical point where superfluidity can no longer 
be quenched by population imbalance, i.e. the CC limit 
becomes 100%. 

In the limit of a BF mixture |3], we observe repul- 
sive interactions between the fermion dimers and un- 
paired fermions. They are parameterized by an effec- 



tive dimer-fermion scattering length of abf = 1.23(3)a. 
This value is consistent with the exact value abf = 1.18a 
which has been predicted over 50 years ago for the three 
fermion problem j4j, but has never been experimentally 
conflrmed. The boson-boson interactions were found to 
be stronger than the mean-fleld prediction in agreement 
with the Lee-Huang- Yang prediction [5 . 

The system for this study is a variable spin mixture 
of the two lowest hyperflne states | t) B.nd \ I) of ^Li 
atoms (corresponding to the \F = 1/2, nip — 1/2) and 
\F = 1/2, mi? = —1/2) states at low magnetic fleld) in 
an optical dipole trap as described in Refs. [HE]. A 
broad Feshbach resonance, located at 834 G [6], strongly 
enhances the interactions between the two spin states. 
The flnal evaporative cooling was performed at 780 G 
by lowering the trap depth. Subsequently, the magnetic- 
bias fleld B is adjusted to a target value with a ramp 
speed of < 0.4 G/ms, changing the interaction strength 
adiabatically. At the end of the preparation, our sample 
was conflned in an effective three-dimensional harmonic 
trap with cylindrical symmetry. The axial (radial) trap 
frequency was lUz/^tt = 22.8 Hz {ujr/l-K = 140 Hz). 

The phase diagram for the fermion mixture was ob- 
tained from the analysis of in situ density proflles of 
the majority (spin 1) and minority (spin J,) components. 
The proflles were recorded using a phase-contrast imag- 
ing technique [2]. Under the local density approxima- 
tion (LDA), low-noise column-density proflles were ob- 
tained by averaging the optical signal along equipoten- 
ital lines (refer to Ref. 2 for a full description of the 
image processing). For typical conditions, the tempera- 
ture of a sample was T/Tpo ^ 0.05, determined from the 
outer region of the cloud [2], where Tpo « 1.0 /iK is the 
Fermi temperature of the majority component measured 
as ksTpo = mti;^_R|/2 (ks is the Boltzmann's constant, 
m is the atom mass, and _R| is the axial radius of the 
majority cloud). 

Figure [T] displays density proflles of imbalanced Fermi 
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FIG. 1: (color online) Density profiles of imbalanced Fermi mixtures in a liarmonic trap. The top row (a-e) sliows tlie averaged 
column density profiles for various magnetic fields (green: majority, blue: minority, black: difference). The black dotted line is 
a zero-temperature Thomas-Fermi distribution fit to the majority wing profile (r > Ri). The middle row (f-j) and the bottom 
row (k-o) show the reconstructed three-dimensional density distributions and the spin polarizations obtained from the profiles 
in the top row. Ri, Ri (dashed dot lines), and Rc (dashed lines) are the radii of the majority cloud, the minority cloud, and 
the superfluid core, respectively. The critical polarizations ctc at the phase boundary r = Rc are indicated by the right arrows. 
The values for _R| (in /xm), Rc/Rf, and Ri/R^ were respectively: for (a,f,k), 381, 0.33, 0.33; for (b,g,l), 380, 0.33, 0.33; for 
(c,h,m), 362, 0.35, 0.59; for (d,i,n), 371, 0.44, 0.72; for (e,j,o), 367, 0.41, 0.76. T/Tfo < 0.05 and Tpo « 1.0 (see the text 
for definitions). 



mixtures for various magnetic fields, showing how the 
spatial structure of a trapped sample evolves in the 
crossover regime. Near resonance, as reported in Ref. [5], 
three distinctive spatial regions are identified: (I) a su- 
perfluid core, (II) an intermediate region of a paritially- 
polarized normal (Npp) phase, and (III) a fully-polarized, 
outer wing. The core radius Rc was determined as the 
peak position in the column density difference profile and 
the majority (minority) radius i?| {Ri) was determined 
from the fit of the outer region, r > i?| (r > Rc) of 
the majority (minority) column density profile to a zero- 
temperature Thomas-Fermi distribution. The local spin 
polarization is defined as cr(r) = (ri| — n|)/(n| -I- n^), 
where rij and are the local majority and minority 
density, respectively. 

Further on the BEC side, the sample has a more com- 
pressed superfluid core, a narrower intermediate normal 
region (gray region in Fig. [T]) and a higher critical spin 
polarization at the phase boundary CTc = cr(i?c). Eventu- 
ally, when B < 780 G, there is no noticeable intermediate 
region, implying that every minority atom pairs up with 
a majority atom in the superfluid core. In Fig. [2] we de- 
termine the critical point for the disappearance of the 
partially polarized normal phase in two different ways. 



Fig. |2ja) shows the phase diagram for the A^pp phase in 
the plane of interaction strength l/kp^a and spin polar- 
ization a. An extrapolation of the critical line to CTc = 1 
yields l/kp^^cO. ~ 0.74(4). Another implication of the 
absence of an iVpp phase is that the size of the minor- 
ity cloud Ri approaches the radius Rc of the superfluid 
core. This extrapolation is conveniently done using the 
dimensionless parameter k = (i?| — i?|)/(i?| — i?^) Hi 
resulting in a value of l/kp^^cO' = 0.71(5). These values 
are in good agreement with recent quantum Monte-Carlo 
(QMC) calculations [SJ. 

The critical point marks the onset of the emergence 
of a BF mixture from the two-component Fermi system. 
One may suspect that near the critical point the equa- 
tion of state of the BF mixture is complex, but we show 
now that a very simple equation of state is sufficient to 
quantitatively account for the observed profiles. Due to 
the external trap potential, the local chemical potential 
varies form zero at the edge of a cloud to a maximum 
value in the center. Therefore, knowledge of the three- 
dimensional density profiles of a single cloud is sufficient 
to obtain the equation of state [51 [31 [HI [IHl [H] ■ 

For a zero temperature mixture of bosonic dimers 
with density rib = ni and mass mb = 2m and un- 
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FIG. 2; (color online) Emergence of a Bose-Fermi mixture 
in the phase diagram for a two component Fermi gas. (a) 
The critical polarization ac as a function of the interaction 
strength l/kpia at the phase boundary. The open circle in- 
dicates the previously measured critical value on resonance, 
aco = 0.36 0. (b) K = (i?^ - Rl)/{R^^ - Rl). The red solid 
lines are (a) an exponential fit and (b) a linear fit to the data 
points. (Jc = 1 and «: = 1 (i.e. i?| = Rc) imply the absence 
of minority fermions in the normal phase. Each data point 
consists of nine to twenty three independent measurements 
and the error bars indicate only the statistical uncertainty. 



FIG. 3: (color online) Characterization of a strongly inter- 
acting Bose-Fermi mixture, (a) The scattering length for 
dimer-fermion interactions abt and (b) for dimer-dimer inter- 
actions Obb in units of the fermion-fermion scattering length 
a. Black solid circles were determined using mean-field the- 
ory and red open circles including the LHY correction for a 
strongly interacting Bose gas. The dashed lines indicate the 
calculated values abf = 1.18a |4] and Obb = 0.6a |13) . Each 
data point represents seven to seventeen measurements and 
the error bars indicate only the statistical uncertainty. 



paired fermions with density rif 



and mass 



TOf = TO, the energy density £ can be decomposed as 
£ = ^-bb + £hi + i^f, where £hh{n\j) and fbf (ub, rif) are 
the boson-boson and boson-fermion interaction energies, 
respectively, and £i = (3/5)cknj^'^ is the kinetic energy 
of fermions [a — (67r^)^/^?i^/2TOf and fi is the Planck's 
constant]. Here we assume that the effective mass of a 
fermion in a dilute mixture is same as its bare mass [12]. 
Under the LDA, the densities n\^{r) and rif (r) in the har- 
monic trap should satisfy 
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where /ifo and /ibo Sire the global chemical potentials of 
fermions and bosons, respectively, referenced to the trap 
bottom. 

For the determination of the boson-fermion scattering 
length abf, we use a mean-field expression for the boson- 
fermion interaction energy £m,m = (27r?i^/TObf )abfn.b7if 
with TObf = TObTOf/(TOb + TOf)=(2/3)m. Since /ifo — 



an 



2/3 



i??/2, Eq. M gives the relation, abf(r) — [^fo(l 



]/(^^^nb). We obtained a value for abf 



by averaging abf(f) over a mixed region (r < Rc with 
nb,f > O.lrio). Here, no is the reference density defined 
as no = (/.ifo/a)^/^. In this analysis, the non-interacting 
outer wing provides absolute density calibration [IV . 

The scattering length ratio Cbf/a turns out to be al- 
most constant over the whole range, 700 G< B <780 G 
where we could study BF mixtures [Fig. |3|a)]. For even 
lower magnetic fields, severe heating occured, probably 
due to molecular relaxation processes. By averaging a to- 
tal of 89 measurements, we obtain Obt = 1.23(3)a, close 



to the exact value Obf = 1.18a calculated for the three- 
fcrmion problem [3] . Our finding excludes the mean-field 
prediction abf — (8/3)a. The detailed behavior above 
750 G requires further investigation. 

We now turn to the determination of the boson-boson 
scattering length abb which parameterizes the boson- 
boson mean-field energy £hh,M = (27r?i^/TOb)abb'^b- ^'^^ 
a given abf , the effective potential for bosons in the pres- 
ence of fermions is Vb(r) = mw^r^ -I- (37r?i^/TO)abfnf (?'). 
Then, Eq. ^ gives /ibo-(27r?^,^/TO)abb7^b(^') = Vb(r). By 
fitting the data in the core region (0.1i?| < r < Rc and 
nb > 0.1 no) to this equation with ^bo and abb as two 
free parameters, we obtained a value for abb- We used 
the value abf determined from the corresponding profiles. 

The effective mean- field values for abb /o. show a strong 
increase by a factor of about 2, as the system approaches 
the critical point [Fig. |3|b)]. We attribute this behavior 
to strong boson-boson interactions causing non-negligible 
quantum depletion in the BEG. In this regime, the equa- 
tion of state has to include beyond-mean-field correc- 
tions, with the leading term given by Lee, Huang, and 
Yang (LHY) as 



LHY — 



2'Kh abb"-b 



2 r 



TOb 



128 J al^n^ 
15 



(3) 



Inclusion of the LHY correction leads to smaller fitted 
values for abb/a, which are now almost constant over 
the whole range of magnetic fields with an average value 
of abb ~ 0.55(l)a. The exact value for weakly bound 
dimers is abb = 0.6a |13j . For kp^ha ~ 1, the LHY cor- 
rection is 0.3fbb,M, i-e. a 30 % correction to the mean 
field approximation. Here, kp^Xi = (67r^nb)^/^. Recently, 
the LHY corrections have been observed via the upshift 
of collective oscillation frequencies for a strongly inter- 
acting BEG [g. 
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FIG. 4: (color online) Observed profiles of strongly interacting 
Bose-Fermi mixtures compared to calculated profiles without 
any adjustable parameter, (a-c) Density profiles of bosonic 
dimers (blue) and unpaired excess fermions (green) for various 
magnetic fields. The numerically obtained density profiles 
(d-f) for bosons and (g-i) for fermions use ttbb ~ 0.6a and 
flbf = 1.18a (dashed: mean-field description, solid: including 
the LHY correction). The horizontal dotted lines in (d-f) 
indicate the boson density corresponding to l/kp.ha — 1. The 
values for Tfo (in mK) and Rc/R-\ were respectively: for (a), 
1.15 and 0.29; for (b), 1.1 and 0.33; for (c), 1.0 and 0.35. 



Our results show that a two-component Fermi mixture 
beyond the critical point can be effectively described as a 
strongly interacting BF mixture. In Fig. |4] we compare 
our experimental data with numerically obtained density 
profiles [15] without any adjustable parameter, showing 
excellent agreement. After including the LHY correction, 
small discrepancies are visible only at the highest boson 
densities exceeding fci?ba = 1, where one would expect 
unitarity corrections. It is surprising that we didn't need 
any beyond-mean-field corrections for the boson-fermion 
interaction. Such corrections have been calculated for a 
system of point bosons and fermions |16j . However, in- 
cluding them into our fit function degraded the quality 
of the fit. Recent QMC simulations have shown that the 
equation of state of a polarized Fermi gas on the BEG 
side is remarkably close to £ — £lhy + £hf,M + £f with 
Obb — 0.6a and Obf = 1.18a down to l/kp^a > 0.5 [8| [T7] 
in agreement with our findings. It appears that the 
beyond-mean-field term is offset by other corrections, 
possibly due to the composite nature of the bosons. Fur- 
ther studies of this rich system could address beyond- 
mean-field terms, characterize the break-down of the BF 
description close to the critical point, and look for finite 
temperature effects. 

One motivation for the realization of BF mixtures is 
to extend studies of '^He-^He mixtures. With tunable in- 
teractions near Feshbach resonances, cold atom systems 
can access a wider regime of the phase diagram. Pre- 
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FIG. 5: (color online) Mean-field phase diagram of a di- 
lute Bose-Fermi mixture. fl{ and jlh are the chemical po- 
tentials for fermions and bosons in units of tt — (6n'^nf)'^^^ 

2mf ^ ' 

and ^^^^^n^, respectively, where fn = and 

nh = ^4^ , "'i'^^s (see Ref. The thick hues indi- 

cate the phase transitions (solid: first-order, dashed: second- 
order). The thin lines represent typical cuts through the phase 
diagram realized in our trapped samples. 



dieted phenomena include phase separation and misci- 
bility [9j [18], boson- mediated, effective fermion-fermion 
coupling [TT, T9 and novel collective excitations PHI [H] • 
The density profiles in Fig. [4] show a sharper boundary 
for higher magnetic fields. This is consistent with Fig. [5] 
which predicts that in the same magnetic field range, the 
transition from full miscibility to phase separation has 
taken place. 

An interacting BF system has been also realized in 
^^Rb-^'OR mixtures [11113]. The ^Li system studied here 
has the advantages of using a single atomic species and 
much longer lifetimes of several seconds, but cannot ac- 
cess attractive boson-fermion interactions. 

In conclusion, a two-component Fermi gas with popu- 
lation imbalance is a realization of a long-lived strongly 
interacting BF mixture. This is a new BF system with 
tunable interactions. Furthermore, it offers intriguing 
possibilities to study the emergence of bosonic behavior 
from a mixture of fermions. 
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